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Abstract. The various schemes for studying rigidly rotating perfect fluids in general 
relativity are reviewed. General conclusions one may draw from these are: (i) There is a 
need to restrict the scope of the possible ansatze, and (ii) the angular behaviour is a valuable 
commodity. This latter observation follows from a large number of analytic models exhibiting 
a NUT-like behaviour. A method of getting around problem (ii) is presented on a simple 
example. To alleviate problem (i) for rigidly rotating perfect fluids, approximation schemes 
based on a series expansion in the angular velocity are suggested. A pioneering work, due 
to Hartle, explores the global properties of matched space-times to quadratic order in the 
angular velocity. 

As a first example of the applications, it is shown that the rigidly rotating incompressible 
fluid cannot be Petrov type D. 
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1 Introduction 

In the past decades, numerous attempts have been made to construct analytic models 
of rotating perfect fluid bodies in general relativity. The static Schwarzschild metric 
of an incompressible fluid ball is known since 1916. The next development, Kerr's 
rotating solution for the vacuum domain, took no less than 46 years to happen. Three 
decades thereafter, the flrst, highly idealized model of a thin rotating disk of dust was 
found by Neugebauer and Meincl[0. There are clear indications from this simple limit 
that the rotating relativistic star model is a task that is not easy to solve. 

A perfect fluid medium in general relativity is characterized by a stress-energy 
tensor of the form 

Tik ^ iii+p) UiUfc - pgik , (1) 

where p is the pressure, p the density and u the four-velocity of the fluid, normalized 
by UiU* = 1. 

The field equations governing the fluid motion are Einstein's gravitational equa- 
tions. 



Rik — \gikR — kTik 



(2) 
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Barotropic fluids additionally obey an equation of state p = p (fi) . The state is 
called stationary and axisymmetric when there exist two commuting Killing vectors, 

it ~ Si — "Sfi ^ 

such that is time-like and space-like. 

The fluid is said to be in circular rotation if the four-velocity u lies in the plane of 
the Killing vectors, 

n = Ni^t + n^^). (4) 

Here N is a normalization factor and is the angular velocity of the fluid. Rigid 
rotation is characterized by the condition that f2 is a constant, and for differential 
rotation: dfl ^ 0. 

Spherical symmetry is characterized by a space-time filling family of 2-surfaces 
dd'^+ sin^-ddLp^ (Delgaty and Lake§, 1998). 

The state is said to possess an equatorial symmetry when it is invariant under the 
substitution i) ^ -k — 'd. 

This paper presents a short review of the subject. The story begins in Sec. 2, with 
a Newman-Penrose presentation of the interior Schwarzschild metric. In Sec. 3, an 
overview is given of the various attempts to create more general models which, then, 
would describe rotation effects. It will be apparent hence, however, that success has 
yet eluded researchers. The main conclusion to be drawn from this situation is that, 
in the interest of better efficiency, there is a need for singling our the few ansatze 
that correctly encapsulate the physics of the rotating system. It is suggested here 
that a very efficient tool for doing this is given by the series expansion in powers of 
the angular velocity fi. This scheme has been worked out by IIartle[|l9| and used 
subsequently for a numerical study of stellar models pOl, 



2 The interior Schwarzschild metric 

The interior Schwarzschild metric, in its pristine form, is written 

ds^ = (A - cosxfdt^ ~ R\dx^ + sin^ x(d6'2 + sin^ 6di^'^)\ , (5) 

where the radius R of the sphere determines the density, /i = and the 

constant A = 3-\/l — r\llV- encodes the radius ri of the star. The pressure of the 
fluid is 

^ 1 3cosx-^ 
^ i?2 A-cosx ' 

The familiar radial coordinate r may be reintroduced by writing sinx = r/R. 

A Newman-Penrose analysis is conveniently carried out by use of a new, advanced, 
time coordinate u such that 

dt^du+ dx- (7) 

A — cos X 
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Then the metric takes the form 

ds^ = {A~ cos xfdu^ + 2i?( A - cos x)dudx (8) 
sin^ x(d6|2 + sin^ OdLp^) . 

This is in a comoving coordinate system since we have 

u^{A-cosx)~^d/du . (9) 

The nuU tetrad is chosen for the metric 

P=— 1 did 

dx ^ R{A~ cos x)du 2R^dx ^ ' 

d id 



\/2i?sinx V^^* smOdip, 

Thus the four-velocity is given by the hnear combination of the tetrad vectors 

Vi^^l + Rn. (11) 

Inserting in the NP equations, we get the Ricci scalar 

A- ^ -4-2C0SX , . 

AR^ A-cosx ^ ' 

and the Ricci-tensor components by 

A _ A s> - ^ ^ 

" A-cosx " ~ 4i?2(A-cosx) ~ ^_cosx ^ ' 

$01 = «'02 = $12 = . 

The spin coefficients satisfy the simple relations r = 7r = a + /? = e — e = 0. 

The Weyl spinor vanishes, ^'if = , where K = 0,1, ...,4 . The interior Schwarzschild 
solution is unique among the perfect fluid metrics in having the property of conformal 
flatness. 

Collinson's theorem (||l,1976) : 

The only conformally flat axistationary perfect fluid is the interior 
Schwarzschild space — time . 



3 An (incomplete) compendium of the ansatze 

In this section we review the various schemes deployed for obtaining stationary and 
axisymmetric perfect fluid solutions of the gravitational equations. The entries are 
arranged in groups characterized by a common ansatz. For a detailed account of the 
state of the subject, cf. 

• Petrov type 
Examples of type D solutions: 

(i) Wahlquist , 1968 solution: This is a rigidly rotating perfect fluid with the equa- 
tion of state /i + 3p = const. 
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(ii) Senovilla]^, 1987 solution: The equation of state is fi — p = const. A later work 
by the same author, Q 1992, presents a wide class of differentially rotating metrics 
containing a free function. 
Relation to type N vacuum: 

Hoenselaers et al.^^, 1994 find a map from the type N vacuum field equations to the 
equations of rigidly rotating perfect fluids. 

• Local Rotational Symmetry (LRS) 

Stewart and Ellis||^, 1968 work out a classification scheme and find some analytic 
solutions. Marklund and Bradley, in a two-part paper [Q, 1997 and 1999, employ the 
curvature representation in their study and obtain examples. 

• Fluid kinematics 

Chinea and Gonzalez- Romero [|l3|, 1992 develop a formalism with differential forms 
adapted to the kinematics of the fluid. 

ChineaQ, 1993 uses this method for getting an analytic solution with differential 
rotation, by assuming that the acceleration of the fluid has a special form. 

• Killing tensor 

The vacuum Kerr metric possesses a nontrivial Killing tensor. This has instigated 
a research of perfect fluid states with a corresponding property by, among others, 
Kramer, Karlovini and Rosquist[^ and Goliath. 

• Vanishing Simon tensor 

Theorem (Kramer [p^, 1985): for rigid rotation, the only possible equation of state is 
of the form + 3p = const. 

Marklund et al.^, 1997 generalize the notion of the Simon tensor for fluids in differ- 
ential rotation. 

• Electric-magnetic Weyl 

Rotation is known to give rise to magnetic curvature. Since the interior Schwarzschild 
solution has zero conformal curvature, it has been surmised that the rotating interior 
Schwarzschild solution will have a purely magnetic curvature. Fodor et al.^^, 1999, 
develop a tetrad formalism in which the electric and magnetic parts of the curva- 
ture explicitly appear. Purely magnetic curvatures are discussed by Lozanowski and 
Mcintosh [|||. 

• Lagrangian symmetry 

Stephani[p^, 1988 investigates the invariance properties of the perfect fluid Lagrangian. 
He finds that the system possesses far fewer symmetries than the vacuum Lagrangian. 

• Static—* stationary symmetry 

Her It pl| , 1988 elaborates on methods by which to a given static state a rotating 
counterpart can be obtained. 

• Geodesic eigenrays 

Lukacs et al[^, 1983, obtain a fluid solution with geodesic eigenrays, exhibiting NUT- 
like behaviour. 

From this impressive array of theoretical efforts, and from the fact that the only 
global solution properly matched to the ambient vacuum region is the inflnitesimally 
thin, flat disk of dust of Neugebauer and Meinel|^, it transpires that the right ansatz 
has not yet been found. 

Hence we may draw two important conclusions for this research: 
CONCLUSION 1: There is a clear need to somehow restrict the scope of 
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possible ansatze. The suggestion here is to resort to perturbation techniques. In the 
next section it will be shown that the approximation scheme developed originally by 
Hartle |19 is an efficient tool for selecting the physically acceptable schemes. 

From the relatively large number of rotating fields with unacceptable causal be- 
haviour, we draw 

CONCLUSION 2: The angular behaviour is a valuable commodity. 



4 Linear perturbations 

In this section we demonstrate the advantages of a perturbative treatment in the slow- 
rotation limit. Although previous works focussed mainly on numerical applications, 
we want to point out that perturbative methods can be used here like a torchlight 
when trying to find the causes of the difficulties met with analytic models. 

(i) The lesson from the vacuum 

Rotating vacuum fields are governed by the Ernst equation |jl^, (1968): 

(ee-l)Ae = 2?v^v^ (14) 

In spheroidal coordinates, this is written 



^ ^x'^ ^l) — - — (I- ^) — 
' dx dy dy 



dx 



2^ 



ix^ - 11 



dx 



Linearizing around the Schwarzschild solution — x, 
^ ^ X + (f 1 real) 

the derivative of the first-order function i^i satisfies the Laplace equation: 



A 



5^6 



dx^ 
Hence 



= 



dx^ 



1=2 



where Pi (y) is the Legendre polynomial of order / and 
^1 = ^ a; J (^j Qi {x) dx^ dxPi {y) + ay . 



(15) 



(16) 



(17) 



(18) 



(19) 
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The analytic solutions of the Ernst equation found by Tomimatsu and Sato ([p^, 
1972) have the rational form 

e = a//3 . (20) 

Fortuitously, the integral ( p^ is of polynomial character. The polynomials in 
have been used for the leading terms in a. 

(a) Perfect fluids 

Hartle (1967) casts the metric of a slowly and rigidly rotating perfect fluid in the 
following form: 



ds^ = e^dt^ - e^dr^ - 



d-d^ + sin"^ ^ (dip - ujdtf + O (n^) . (21) 



From the large-r behaviour it follows that the rotational potential w is a function of 
the radius alone, lu — uj (r). 

From the Schwarzschild solution (^ for an incompressible fluid we get: 

e'' = (A-cosx)^ e^ = -^. (22) 

sm X 

The potential lo satisfies Heun's equation 

{A - z){z^ _ 1)^ - (3^2 _ 5Az + 2)4^ +iAuj = (23) 
dz^ dz 

with z = cos X ■ 

This equation has proved to be a stumbling block in the way of progress for three 
decades. Naive attempts at removal do fail. For example, one may ask for the effect 
of coordinate change. Since the uj equation is a special case of the six-parameter Heun 
Eq., it would seem conceivable that in a better coordinate system, solution is more 
amenable. Such a new coordinate system is provided, for example, by the property 
of the interior Schwarzschild metric that it is conformally flat. As a consequence, the 
metric can be conformally mapped to the Einstein Universe (as the case is with the 
Minkowski spacetime): 

ds'^ = n'^ [dt^ - dx^ - s\r? x{di}^ + sin^ i9d(p'^)] . (24) 

The conformal factor is il^ = R^{A^ — 1)/{A — cosx)'^. The transformation to the 



conformal form (24) is as follows, sina; = {A"^ — iy^^{A — cosx)^^sinx ■ We find, 
however, that the perturbation equation is of the Heun type again. 

In a like manner, attempts to change the dependent variable have failed to turn 
up a solvable equation. A straightforward thought here is to try and use the Ernst 
potential in place of the rotation potential lu in the metric 

ds'^ = f{dt + ujdipf (25) 
-r' [e'^dr' + (1 - r^) ((1 - y^'dy' + (l - y') V)] , 
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where y — cos-d. The field equation for the Ernst potential ^ is 

{a - 1) AC - 2CVe • VC = ifc (CC - 1) + 1)' (m + 3p) (26) 

with p the pressure and the density of the fluid. 
Integration of the Bianchi identity yields 

(m + P) f^'^ = const. (27) 
For the interior Schwarzschild metric, 7 = 0, and 

3r — /I 

f = iA-r)\ fcM = 3, kp=- . (28) 

A — r 

Let the rotational perturbation be introduced, 

e-^+^Ci- (29) 

The linearized field equation for is separable; the y dependence satisfies the Leg- 
endre equation of index I . 

It is possible to integrate the radial equation for / = 0: 

1 ^ \r'^+3Ar + A^r-2-3A'^ / r-^ 

6 = oiCi-Ca 7^= +3Aln r+V^-^-l }. 

This would appear quite encouraging. However, from asymptotic matching, one finds 
that ^1 is proportional to y, and hence we need the integral for I — I, which is not 
known in a closed form. 



5 Two theorems 

In this section we describe how to restrict the scope of choice in order to meet Con- 
clusion 1, by applying the perturbative treatment to various systems. 

5.1 Theorem on incompressible Quids 

To second order in il, the metric of the rigidly rotationg incompressible fiuid is|pj| 

ds^ = {l + 2h){A-cosxydf ~ R^{l + 2k)sin'^Xsin^Hd(p-ujdtf 

-i?2(l + 2m)dx^ - i?2(l + 2n) sin^ xrfi9^ . (30) 

Here lo is of first order, and h, k, m, n of second order in the angular velocity Q. 
A comoving tetrad may be chosen 

f 1 — h , ^, 9 i?2 sin^ iJsin^ Y 
4 -; + - ^ 

A-co^X 2(A-cosx) 

Gn = —. , — (31) 

" cosx ^ R ^ ' 

^ 1 — n ( _Rsin??sinx(ti; — i7) 



-Rsinx ' (A - cosx) 

- Vl 



V _ k-l ^ f 2 o2^ -Rsini^sinx 



i? sin 7? sin X 2 (A -cosx) 



2 
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The vorticity is given by the first-order rotation coefficients 



1 



UJ2 



2 (A- COS x) 
sinz? 

2 (A- COS x)' 



dui 

2 cosd{Q, — bj) — sini9— — 
ov 



2{ijj — il)(ylcosx — 1) + sinx {A — cosx) 



duo 
dx 



(32) 
(33) 



The changes in the other rotation coefficients are of the second order. 
The magnetic part of the Weyl tensor is of order fi, 



Hi 

H3 



cos-d 



du) 



R (cos X — A) dx 



Ho 



1 



-Hi 



(34) 



^2 COS X ^ sin ■& Lo — ^ 



R sin X R sin x A — cos x 

The eiectric part Ei is of order fi^. The condition for the space-time to be Petrov 
type II, to quadratic order, is H^ = and Hi = i?2- Hence the Petrov type of this 
incompressible fluid cannot be II. 

The Petrov type is D if i?| - 2Hl - 2i7| - Si/ii/a = 0. The general solution of 
this condition has the form 



Lo = C(.4/cosx- 1)^ 



(35) 



This is not a solution of the Heun equation (g3|) which proves the following 

THEOREM [p^ : A rigidly rotating incompressible fluid ball in an asymptoti- 
cally flat vacuum exterior cannot be Petrov type D. 



5.2 Theorem on the Wahlquist sokition 

The Wahlquist metric is an example of type D solutions. It is given by| 



ds^ = 



f {dt - AdifY - {e+e) 



{i~k^e)hi 



hih2 



-dif^ 



where 



/ 

hi (C) 
h2 {£,) 



hi 



1 + + 4 



A = ro 



ehi+eh2 ,2 

hi - h2 



C-- VT^kV arcsin (kC) 



C - f Vl + k^^arcsinh (k^) 



The constant is defined by h2 {^a) = . 
The pressure and density are 

p = i^io (1 - K^f) M = 5M0 (3kV - 1) 



(36) 



(37) 
(38) 
(39) 

(40) 
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The equation of state is 
2k2 



(41) 



THEOREM [27 1 : By quadratic approximation in , 

The Wahlquist metric cannot be matched to an asymptotically flat vacuum exterior. 



6 How to secure the angular behaviour? 

I do not reaUy have a satisfactory answer to Conclusion 2 concerning the angular be- 
haviour. Ah I can offer here is a particular recipe by which one creates a relativistically 
rotating version of a certain static system. There is no guarantee, however, that the 
space-time so obtained will have the desired type of matter source. Succinctly stated, 
the idea is as follows. The constant-time sections of the Schwarzschild fluid ball (|^) 
are 5*"^ spheres. Now an may be squeezed and simultaneously set in rotation about 
the z axis following these steps: 
— Embed in E"^ with the metric 

ds^ = {dx^f + {dx^f + {dx'^f + {dx^f . (42) 

Introduce polar coordinates (i?, x, 

x^ + ix^ = i? sin X sin de^"^ 

x^ = R sin X cost} (43) 
x'^ = R cos X ■ 

The surface is given by i? = const. 

— Squeeze and set simultaneously in rotation by putting 

X-^X + ia (44) 

in x^ + ix^ . Thus the R = const, surface is flattened, and has the metric 

ds^ = (sin^ X + sinh^ a cos^ d){d'd'^ + sin^ Mip"^ ) 
-|- sinh a cosh a sin^ dxdf — sinh'^ a sin^ ddtp^ . 

— Embed in the Einstein Universe. For R—\, the metric is of the form 

ds^ = dt^ - dx^ - sin^ xid^"^ + sin^ M^p^) . (45) 

— Soak up the slack dx dip term by introducing the advanced time coordinate 

u^t + x- (46) 
The space-time metric then takes the form: 

ds^ = (dii — sinh a sin^ ?9diy9) ^ 

— 2 (dii — sinh a sin^ i}d(p^ (cosh adx ~ sinh a sin^ "ddip) 
- (sin^ X + sinh^ a cos^ -d) {dd'^ + sin^ 'ddip^ ) . 
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This metric represents a squeezed and rotated Einstein universe. In the Hmit as 
i? — > oo, we recover Minkowski space-time. 

CONJECTURE : By continuity in the rotation parameter a, 
There is an open neighbourhood of Minkowski space-time in the minisuperspace of 
squeezed Einstein universes that is causally well-behaved. 

Apphcation of this procedure to the interior Schwarzschild metric yields rotating 
models with various stress-energy tensors, depending on details of the embedding. 

It is a pleasure to thank the Organizers of the Journees Relativistes Conference for providing 
a seamless platform for scientific activity and for hospitality. Part of this work has been done 
under the contract OTKA T022533. 
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